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in ■ ABSTRACT 

O . 

^ I We study the deconfinement phase transition and monopole prop- 

I erties in the finite temperature 3D compact Abehan gauge model on 

^ ' the lattice. We predict the critical coupling as function of the lattice 

p • size in a simplified model to describe monopole binding. We demon- 

^ ■ strate numerically that the monopoles are sensitive to the transition. 

• • , In the deconfinement phase the monopoles appear in the form of a 

• ^ ! dilute gas of magnetic dipoles. In the confinement phase both mono- 

^ I pole density and string tension differ from semiclassical estimates if 

I monopole binding is neglected. However, the analysis of the mono- 

pole clusters shows that the relation between the string tension and 
the density of monopoles in charged clusters is in reasonable agree- 
ment with those predictions. We study the cluster structure of the 
vacuum in both phases of the model. 



1 Introduction 

Compact Abelian gauge theory in three Euclidean dimensions is a case where permanent 
confinement is proven and qualitatively understood 0. In order to gain some ex- 
perience for more realistic theories, it is interesting to study how this mechanism ceases 
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to work under special conditions. High temperature is such a case. In this paper we 
are going to revisit the finite temperature deconfining phase transition. We will empha- 
size the aspect of monopole binding which explains the breakdown of confinement. In a 
companion paper, we extend these studies to the case of non-vanishing external fields. 

Compact QED theory possesses Abelian monopoles as topological defects appearing 
due to the compactness of the gauge group. Considering the 3D theory as the static 
limit of a 4D theory, the monopoles are just magnetic monopoles at rest, and the com- 
ponents of field strength are magnetic. In a three dimensional theory, the monopoles 
are instanton-like objects: instead of tracing world lines (as they do in AD) they oc- 
cupy points. The plasma of monopoles and anti-monopoles can explain the permanent 
confinement of oppositely charged electric test charges [|l| in bound states, kept together 
by a linear potential. In the language of magnetostatics, confinement appears due to 
the screening of the magnetic field induced by the electric current circulating along the 
Wilson loop. Monopoles and anti-monopoles form a polarized sheet of finite thickness 
("string") along the minimal surface spanned by the Wilson loop. The formation of the 
string (observed in the lattice simulations in Ref. [Q) leads, for non-vanishing electric 
current, to an excess of the free energy proportional to the area. 

At finite temperature the phase structure becomes non-trivial. What we have in 
mind, is compactifying 3D (2 + 1)D in the "temporal" (third or z) direction. The con- 
finement-deconfinement phase transition was studied on the lattice both analytically 0] 
and numerically [0. According to the Svetitsky-Yaffe universality arguments P, an 
interpretation of the transition has been attempted in terms of the U{1) vortex dynamics 
of the corresponding 2D spin system. The phase transition — which is expected to be of i 
Kosterlitz-Thouless type — was demonstrated to be accompanied by restructuring of 
the vortex system 0. The vortices are described by a 2D U{1) spin model representing 
the dynamics of the Polyakov line (see also the discussion in Ref. [Q). Approaching the 
transition temperature, vortices and anti-vortices start to form bound states. In the high 
temperature phase no unbound vortices and anti-vortices are left. 

In the present paper we discuss an interpretation of deconfinement starting from the 
confinement picture outlined above, in terms of magnetic monopoles. The confining 
plasma of the monopoles and anti-monopoles turns into the dipole plasma at the de- 
confinement phase transition. The dipole plasma is inefficient to completely screen the 
field created by the electric currents running along the pair of Polyakov lines. In this 
case the screening mass vanishes while the magnetic susceptibility of the "medium" is 
smaller than unity. Both monopole and vortex binding mechanisms of the deconfinement 
phase transition have been discussed for 3D finite temperature Georgi-Glashow model in 
Refs.|g] and [^], respectively. 

In the finite temperature case, strictly speaking, there is a problem to call all the fields 
"magnetic" as we did above when we summarized the zero temperature case. Similar to 
Ref. 0, the confinement aspect itself will be illuminated in terms of the U{1) valued 
Polyakov lines in the third direction and of Polyakov line correlators representing pairs of 
charges separated in 2D space. In {2 + 1)D there is no symmetry anymore between the 
three components of the field strength tensor. The closest relative of the true magnetic 
field is Fi2 distinct from the others, while there is still a symmetry between F13 and 
F23. With this distinction in mind one can conditionally call them the "magnetic" and 
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"electric" components of the field strength tensor, respectively. As long as one does 
not introduce external fields, even at finite temperature there is no need to distinguish 
between them. The sources of the respective fluxes will be simply called "monopoles" or 
"magnetic charges" in the following. 

The binding of the monopoles is not isotropic. It happens mainly in the 2D space 
direction due to the logarithmic potential between the monopoles separated by a large 
spatial distance. As a consequence of the periodic boundary conditions in the temporal 
direction the force between the monopoles and ant i- monopoles vanishes at half of the 
temporal extent. Therefore the potential in the temporal direction is weaker than in 
the spatial directions. As a consequence, the spatial size of the monopole bound state 
is expected to be smaller than the size in the temporal direction. This means that the 
dipoles are dominantly oriented parallel or antiparallel to the S"^*^ direction. The monopole 
deconfinement scenario raises the question whether the monopole properties (such as 
pairing and orientation) could be influenced by an eventual external field. This aspect 
will be addressed in a companion paper. 

With or without external field, the deconfining mechanism by monopole pairing seems 
to have interesting counterparts in more realistic gauge theories. The formation of mono- 
pole pairs is qualitatively similar to the binding of instantons in instanton molecules with 
increasing temperature in QCD suggested to be responsible for chiral symmetry restora- 
tion [|lOl. In the electroweak theory, the formation of Nambu monopole-anti-monopole 



pairs, a remnant from a dense medium of disordered Z-vortices and Nambu monopoles 
which characterizes the high-temperature phase, is accompanying the transition towards 



the low-temperature phase ||TT|. Note also, that a dipole vacuum, although not confining. 



still has a non-perturbative nature [|T2 



The plan of the paper is as follows. In Section |^ we estimate the critical coupling of 
the confinement-deconfinement phase transition based on a monopole binding model for 
a finite lattice. In the next Section the transition is numerically located for a lattice size 
32^ X 8 and confinement properties are studied. We present various monopole properties 
including dipole formation based on a cluster analysis in Section ^. We study the relation 
of the monopoles and dipoles to the phase transition in Section |^. We briefly summarize 
our results in the last Section. 



2 Some heuristic considerations 

In ?)D compact electrodynamics there are monopoles interacting via the Coulomb poten- 
tials, 

S =^^qaqbVT{Xa- Xb) , (1) 

a, 6 

where qa and Xa are, respectively, the charge (in units of the elementary monopole charge, 
g„i = 27r/(yf3, where is the three dimensional coupling constant) and the 3D position 
vector of the a^^ monopole. The subscript T indicates that the interaction potential Vr 
eventually depends on the temperature. 
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At zero temperature the monopoles are randomly located in the Euclidean TZ^ space 
and the classical interaction potential between the monopole and anti-monopole is in- 
versely proportional to the distance R between the objects, Vo(-R) = — (47ri?)~^. At finite 
temperature T the monopoles live in the TZ^ x Si space (with Si being a circle of perimeter 
T~^) and the interaction potential gets modified. At small separations between monopole 
and anti-monopole the interaction is zero-temperature like, Vt(x, z) = Vo{\/x^ + z^) + . . ., 
where x = {x,y,z) = (x, z). At large spatial separations x the potential between 
monopoles is essentially two-dimensional 

1^t(x,2) = -2Tln|x| + ... , |x|T>l. (2) 

However, the interaction between monopoles separated by a distance z in the third (tem- 
perature) direction is of the 3D Coulomb type for small spatial inter-monopole distances, 
|x|T -C 1: Vt(x, z) = — (47rz)~^, zT -C 1. The force between monopoles and anti- 
monopoles at a distance z = 1/ (2T) vanishes due to periodicity in the temperature direc- 
tion. Thus one might expect that at finite temperature the monopoles form magnetically 
neutral states which are bounded in the spatial directions. However, the dynamics of the 
monopoles in the temperature direction is not restricted by a logarithmic potential. 

Thus, at zero temperature the system exists in the form of a Coulomb gas of magnetic 
monopoles and anti-monopoles. In this phase the medium confines electric charges 
As temperature increases, the three-dimensional Coulombic potential turns into a two- 
dimensional logarithmic potential for spatial monopole interactions. The monopoles and 
anti-monopoles become weakly confined and form more and more dipole bound states. 
The dipoles have a finite average spatial size (the distance between the magnetically 
oppositely charged constituents) which is a decreasing function of the temperature since 
the interaction potential between the particles rises as temperature increases, cf. eq. (H). 

In the low temperature regime, this dipole size would be still larger than the average 
distance between the particles inside the plasma, and therefore only a small fraction 
of monopoles residing in actual dipoles is mixed with an weakly correlated monopole- 
anti-monopole component. At sufficiently large temperature, however, the typical dipole 
size becomes smaller than the interparticle distance in the plasma and the system turns 
into a pure dipole plasma. The confinement property is closely related to the Debye 
mass generation effect which is absent in the pure dipole plasma [|T^. As a consequence, 
the confinement of electrically charged particles disappears. The system experiences a 
confinement -deconfinement phase transition due to the monopole binding mechanism. 

One can use these heuristic arguments to estimate the phase transition temperature. 
In continuum theory this analysis was done in Ref. |^ where compact electrodynamics 
was represented as a limit of the Georgi-Glashow model. The phase transition in this 
theory happens at a temperature T = g^/{27[). This result has been obtained under the 
condition that the average size of the effective magnetic dipole is not an infrared divergent 
quantity as it is the case in the confinement phase. 

However, in lattice gauge theory the considered quantities are all finite and the con- 
siderations should be modified compared to the continuum case. The difference between 
monopole and dipole plasmas can only be seen if the mean distance f between the con- 
stituent monopoles becomes comparable to the dipole size d. The distance f can be 
expressed via the density of the monopoles p, as f = p~^/^. Thus, the phase transition 
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happens when the dipole size and the average distance between monopoles become of the 
same order, 



d = ip 



-1/3 



(3) 



where ^ is a geometrical factor of order unity. For both quantities, d and p, estimates can 
be easily obtained on the lattice while the factor ^ is to be defined from a simulation. 

We consider the 3D compact U{1) gauge model on the Li x Lt lattice with the action 
written in the Villain representation: 



Z=jve exp{-^||d^ + 27™||^}, 



(4) 



n(c2)e^ 



where 9 is the compact U{1) gauge field and n is the integer-valued auxiliary tensor field 
variable. (3v is the Villain coupling constant. 

To relate this to the numerical simulations, we also consider the formulation of the 
compact U{1) gauge theory with Wilson action: 



S 



cos 6 1 



(5) 



The Villain coupling constant (3v is related to the Wilson coupling (3 as follows 

-1 



2 log I 



(6) 



where /o,i are the standard modified Bessel functions. 

The partition function (H) can be rewritten in the following "grand canonical" form, 
i.e. represented as a sum over monopole charges in the (dual) lattice cubes 



Z oc Zmon = E exp|-27r^/3y(m, A ^m)} . 

r?i{c3)6^ 



(7) 



Here is the inverse of the Laplacian operator on an asymmetric lattice, rric denotes 
the monopole charge in the cube C3. The inverse Laplacian for lattice sizes Lg, Lt is given 
as follows: 



A ^(x]Ls,Lt) 



2L2L,^5^3-EiiCosp,' 



where pi,2 = 0, . . . , 27r(Ls - and = 0, . . . , 27r{Lt - 1)/ Lt. 

In order to estimate the average distance between the monopole and anti-monopole 
constituents in a dipole state {i.e., the dipole size) we use the "canonical" monopole-anti- 
monopole (dipole) partition function which can be easily read off from eq. (|^: 



4^^ = const - J2 exp{An^M^-\x;L,,Lt)-A-\0;L,,Lt)]} 



(9) 
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the sum extents over all lattice separations x between monopole and anti-monopole. The 
zero distance between these objects is excluded (since this case does not correspond to a 
dipole state). The r.m.s. dipole size J is given by 

(?(/3y; L„ Lt) = ^^Y. exp{47r2/5y [A~\x- L„ Lt) - A-i(0; L„ U)] } , (10) 

where the actual distance squared, x^, is evaluated taking into account the periodic bound- 
ary conditions of the lattice. The sums cannot be taken analytically. 
The monopole density p can be read off from eq.(^, 

p(/?;L„LO = 2exp{-27r2/3v'(/?)A-i(0;L„L0}, (11) 

where the dependence on the lattice geometry is indicated explicitly. Note that in this 
formula no interaction between monopoles is taken into account and we refer to it as to 
"bindingless" . Only the local "self-interaction" of monopoles is accounted for via the 
Coulomb propagator A~^(0) in the fugacity. We are discussing the binding effects on the 
monopole density in Section ^. 

The geometrical factor ^ is to be defined from the numerical data. To this end we 
assume that this factor is a constant i quantity which does not depend on the lattice 
extensions. Indeed, it gives an estimate how large the intra-dipole distances should be 
compared to the monopole density in order to have the dipole field screened. This is a 
quite strong assumption which, however, turns out to be reasonable, as will we see below. 
To define the factor ^ we substitute eqs. (|10D and (|ll|) into eq. (|) and use numerical values 
for [3c presented in Ref. 0. For the lattices 16^ x Lt, Lt = 4,6,8 we get, respectively: 
^ = 0.723(58), 0.622(47) and 0.646(116). These numbers coincide with each other within 
numerical errors. Taking the average over Lt we get ^ ^2/3. In what follows we take 

e = 2/3, (12) 

and then solve eqs.(^|TT],|IU]) with respect to the Villain coupling Pv Then we finally 
estimate the critical Wilson coupling with the help of eq. (|^). The results for lattices 
of various sizes are represented in Table |I| and compared with pseudocritical couplings 





Ls 


= 16 


Ls 


= 32 


Ls = Q^ 


Lt 


f3f 


/3c 




/3c 


/3f 


4 


1.87 


1.83(2) 


2.01 




2.10 


6 


2.04 


2.08(2) 


2.26 


2.18(3) 


2.44 


8 


2.12 


2.14(5) 


2.39 


2.30(2) 


2.62 



Table 1: The critical coupling constant /?*'* calculated using eqs. (|6|,|3| JT0|JTT|JT^) for different 
lattices Li x Lt compared to lattice Monte Carlo results of Ref. 0. Note that our results 
for the lattice 32^ x 8 are slightly higher than that of Ref. [0, see forthcoming Sections. 

(3c obtained in lattice simulations of Ref. 0. The agreement between the data and our 
estimates is within 4%. Thus the simple heuristic arguments based on the monopole- 
dipole picture work surprisingly well. 
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3 Phase transition and confinement 



We have performed our numerical study of (2 + 1)D compact electrodynamics using the 
Wilson action (^). The lattice coupling /3 is related to the lattice spacing a and the 
continuum coupling constant gs of the 3D theory as follows: 

(3=^, (13) 

a 93 

Note that in three dimensional gauge theory the coupling constant gs has dimension 
mass"^/^. 

The lattice corresponding to the finite temperature is asymmetric, x Lt, Lt < 
Lg. In the limit — > oo the "temporal" extension of the lattice Lt is related to the 
physical temperature, Lt = l/(Ta). Using eq. (p!3| ) the temperature is given via the 
lattice parameters as follows: 



L- £. 

93 Lt 



(14) 



Thus, at fixed lattice size lower (higher) values of the lattice coupling constant P corre- 
spond to lower (higher) temperatures. 

Our simulations have been performed mainly on a 32^ x 8 lattice. We do not intend 
to study in the present paper finite size scaling aspects of this model. The local Monte 
Carlo algorithm is based on a 5-hit Metropolis update sweep followed by a microcanonical 
sweep. For better ergodicity, in particular in the presence of an external field (considered 
in a companion paper), also global updates are included. Following the ideas of Ref. il5| , 
the global refreshment step consists in an attempt to add an additional unit of flux with 
randomly chosen sign in a direction randomly selected among the three, to the dynamical 
gauge field subject to a global Metropolis acceptance check. 

For example, one unit of flux in ij plane is introduced with the help of the following 
gauge field shift ^ ^ {^^ + ^^Lod2n■ 



~ TT 

6j = —{2xi — Li — 1) , Oj = for xj ^ Lj , 
Li 

~ 2ti 

Oi = ——{l-Xj), ek = 0, ky^i,j. 

Li Lj 

The acceptance rate of the global step changes within the considered (3 range from roughly 
0.7 (confinement phase) to 0.2 (deconfinement phase). One total Monte Carlo update 
cycle consists of two combined local Metropolis and microcanonical sweeps (requiring an 
acceptance rate of 0.5 for the Metropolis step) and the global update described above. 

In order to localize the deconfinement transition, it is convenient to study the expec- 
tation values of the two bulk operators, 

(l^l) = 4(IE^WI)' (l^r) = 4(IE^WI')' (15) 

constructed from the Polyakov loop, 

Lt 

L(x)=exp{z^^3(x,z)}, (16) 
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here x = {x,y) is a two-dimensional vector. In the deconfinement phase the quantity \L\ 
is of the order of unity, while in the confinement phase it is close to zero in a finite volume 
and vanishes in the infinite volume limit. 

The behaviour of the expectation value of the Polyakov loop vs. lattice coupling (3 is 
shown in Figure |l](a). The low temperature phase, j3 < (3c corresponds to the confinement 




Figure 1: (a) The expectation value of the absolute value of the average Polyakov loop ([T5|) 
and (b) its susceptibility (O) vs. [3. 



phase, while the high temperature phase is deconfining. 
The susceptibility of the Polyakov loop 

XL = (1^1')- (1^1)' (17) 

is shown in Figure 0(b). The peak of the Polyakov loop susceptibility corresponds to the 
pseudocritical (3c of the deconfinement phase transition. We have fitted the susceptibility 
near its maximum by the following function: 

^2 



x!H/5) 



'1^ 



ci + (/3-/3c)2' 

where the critical coupling was estimated to be (3c = 2.346(2) which is quite close to the 
result of Ref . [0 . The best fit is shown in Figure |l|(b) by a solid line. 

To calculate the string tension we use "plane-plane" correlators of two Polyakov loops. 
In addition, averages of temporal Wilson loops have been studied, too. More precisely, 
in (2 + 1)D, we define first sums of the Polyakov loops along a line parallel to a spatial 
lattice axis {e.g. in the y direction): 

Ls 

Lpi^neix) = Y.Hx^y) ■ (19) 
y=l 

The correlator of the plane-plane correlators may be written as a sum of point-point 
correlation functions, 

(I^piane(0)L;i_(a;)) = I: (L(0,yi)L+(x,y2)). (20) 

2/1,2=1 
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The form of this correlator is expected to be: 

(^piane(O)Lpi^ne(x)) = const ■ cosh aLt (^^ " y 
where a is the "temporal" string tension. In Figure ^(a) we show the fit of the Polyakov 
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Figure 2: (a) Fit of the Polyakov plane-plane correlator ( pUD using eq. ( PT| ) in the con- 
finement, (3 = 1.8, and deconfinement, (3 = 2.5, phases, (b) String tensions as functions 
of P compared with the bindingless theoretical result (^). 

plane-plane correlator (|20| ) by this fitting function in the confinement {(3 = 1.8) and 
deconfinement (/5 = 2.5) phases, respectively. 

In Figure ^(b) we present the fitted string tensions as function of (3. Above (3c the 
string tension quickly drops down but stays non-zero due to finite volume effects. The 
temporal string tensions obtained using either the Polyakov loop plane-plane correlators 
or the temporal Wilson loop averages roughly coincide with each other. The dashed curve 
represents the theoretical prediction for the string tension |T^, |TB[: 



:22) 



Agreement between the prediction and the numerical results is reached only in the 
vicinity of the phase transition point, (3 ~ 2.3. In order to understand these differences, 
we turn now to a closer investigation of the monopole properties. 



4 Properties of the monopole— anti— monopole system 

The basic quantity describing the behaviour of the monopoles is the monopole density, 
P = J2c\^c\/{L^Lt), where nic is the integer valued monopole charge inside the cube c 



defined in the standard way [jT7| 



27r 



El 

P&dc 



-If [6 



-PJmod27r ' 



(23) 



where the plaquette orientations relative to the boundary of the cube are taken into 
account. The density of the total number of monopoles is a decreasing function of the 
lattice coupling (3 (or the temperature) as it is shown in Figure ^ by circles. At high 
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Figure 3: The density of all monopoles and of monopoles in neutral clusters vs. (5. 



temperatures (large (3) the monopoles are dilute and form dipole bound states. Typical 
monopole configurations in both phases are visualized in Figure ^. 




(a) (b) 

Figure 4: Typical monopole configurations for (a) the confinement phase {j3 = 1.6) and 
(b) the deconfinement phase (/3 = 2.5) . 



In Figure || we show by the dashed line the density of the monopoles calculated using 



eq. ([TT|) for comparison. As in the case of the string tension, the predicted monopole 
density is in agreement with the numerical data only near (3^- 

Equation (0) is based on the single monopole contribution to the partition function, 
thus it does not take into account pairing of the monopoles. The effect of the constituent 
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monopole pairing (dipole formation) due to finite temperature can explain the deviations 
from the bindingless case seen in this Figure. In the confinement phase the density of the 
monopoles is larger than the prediction of eq. ([TT|) . Indeed, we expect that the formation 
of the bound state decreases the total energy (action) of the chosen monopole and anti- 
monopole. As a result binding favours the creation of additional monopoles by quantum 
fluctuations. This tendency increases with larger (3, however the cost of creating new 
monopoles grows, too. 

Note that the entropy of the bound state is smaller than the entropy of a free monopole 
and an anti-monopole. However the entropy effect does not seem to change essentially 
near the phase transition. 

We remind the reader that on the classical level the dipoles are formed both in the 
confinement phase and in the deconfinement phase due to logarithmic potential between 
the monopoles. However, at low temperatures the dipole size is larger than the average 
distance between the monopoles and, therefore, the dipole formation does not destroy 
confinement. 

One can analyse the monopole pairing studying the cluster structure of the monopole 
ensemble extracted from the Monte-Carlo configurations. For our purposes, clusters are 
defined as follows: clusters are connected groups of monopoles and anti-monopoles, where 
each object is separated from at least one neighbour belonging to the same cluster by a 
distance less or equal than -Rmax- In the following we use -R^ax = 3 which means 
that neighbouring monopole cubes should share at least one single corneiQ. Note that 
the increase of the coupling constant leads not only to an increase of the temperature, 
eq. (0), but to a decrease of the lattice spacing a as well, eq. (p!3[) . Thus at different j3 
the same characteristic distance -Rmax corresponds to different physical scales. Therefore 
our results below are of only quahtative nature. 

A monopole cluster is neutral if the charges of the constituent monopoles sum up to 
zero. We show the density of monopoles belonging to neutral clusters as a function of (3 in 
Figure ^ by squares. The difference between this density and the total monopole density 
amounts to a factor three at /5 ~ 1 and becomes smaller at larger (3. At large (3 (entering 
the deconfinement phase) approximately every second monopole or anti-monopole belongs 
to a neutral cluster. At still larger /3's almost all monopoles are in neutral clusters. 

We are confident that the fluctuation of monopole numbers signal the deconfining 
phase transition. This is demonstrated studying the second and (modified) fourth Binder 
cumulants of the total number of monopoles and anti-monopoles [and of the number of 
(anti-)monopoles being part of neutral clusters]. We present in Figure ^ the cumulants, 
with M denoting the respective number : 

R - ^^'^ 1 (24) 

= 25 

Similarly to the Polyakov line susceptibility these quantities are suitable to localize the 



^ In Rcf . a similar definition has been used to investigate tightly packed clusters with i?„ 
In our case the condition for the cluster is more relaxed. 
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Figure 5: The second (a) and the fourth Binder (b) cumulants according to ( p5|) for the 
total (anti-)monopole density and the corresponding densities enclosed in neutral clusters. 
The fits are shown as solid lines. 

deconfining phase transition. We fit these cumulants by 



n = 2,4 



(26) 



The fits are shown by solid lines in Figure ^ and the results for the pseudocritical couplings 
(3c are given in Table ^ 



cumulant 


2nd 


4th 


total 


2.380(3) 


2.404(4) 


neutral 


2.379(5) 


2.372(3) 



Table 2: Pseudocritical couplings jSc from the fits to the Binder cumulants (|25| , p6D 



Some details on the cluster structure at various values of P can be seen in Figure |^(a). 
We show the fraction of monopoles and anti-monopoles being part of clusters of size 
A^. The cluster size is the number of monopoles and anti-monopoles which belong to the 
given cluster. There is no separation according to the cluster's charge. In the confinement 
phase, (3 = 1.5, the fraction of monopoles is slowly decreasing with the cluster size A^. 
The percentage of isolated (anti-) monopoles (A^ = 1 clusters) amounts to roughly 45 % 
while clusters (with a size up to = 10) contain the rest. 

At the phase transition point [jS ~ 2.3) the number of (anti-)monopoles enclosed in 
larger clusters drops drastically. The monopole vacuum is composed mostly of individual 
(anti-)monopoles (60 %) and dipoles (40 %). This observation can be reconciled with our 
theoretical expectation that all monopoles must become paired only if we accept that the 
"unpaired" monopoles are actually part of dipoles of size bigger than -Rmax- Deeper in 
the deconfined phase, however, at /3 = 2.8 practically 90 % of the (anti-) monopoles form 
tight bound states with sizes smaller than -Rmax = ci. 
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Figure 6: (a) The cluster structure at various of the couphng constant (3. Cluster distri- 
bution is shows as a function of the number of constituent monopoles inside clusters, A^. 
(b) The cluster shape function, eq. ([27|), for various [3. 



As we have discussed above, we expect that the force in the spatial directions is larger 
than the force along the temporal direction z. This fact can be qualitatively analysed 
with the help of the following "cluster sphericity" : 

R.UN) = I ^'f , , (27) 

where (|Ax|)^ is the average distance from the center of the cluster in x-direction etc. 
for cluster size N . If the clusters are elongated predominantly in the temporal direction 
this quantity would be larger than unity, and smaller otherwise. We show the dependence 
of the sphericity i?size on the cluster size for various /3 values in Figure |§(b). Small 
clusters are directed predominantly along the temporal direction, as expected, at all [3. 
With larger (3 the elongation becomes stronger. For large clusters the direction of the 
cluster is random, since in this case the cluster shape function is very close to 1/ a/2 (this 



directly follows from the definition (|27|) ). This random limit is marked by the solid line 
in Figure |. 



5 Confinement and monopoles 

We have observed that agreement between predictions from a theory without monopole 
binding and the finite temperature simulation results is reached only in the vicinity of 
the phase transition point, j3 ~ 2.3. In the confinement phase both measured temporal 
string tension and monopole density are larger compared to the bindingless predictions, 
see Figures H(b), ||. 

As we have discussed, it is due to monopole binding that the density of the monopoles 
is increased compared to the non-interacting case. However, the size of the dipoles in the 
confinement phase is larger than the average distance between the ordinary monopoles 
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calculated from their total density. Therefore, the monopoles bound in dipoles due to 
classical logarithmic potential still give a contribution to the string tension. 

It is interesting to check how the monopole density fits into the theoretical predictions 
of the string tension (0). Using that predicted relation, we compare in Figure ^ the ratio 
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Figure 7: The ratio of the temporal string tensions (|28|) vs. (3. 

R between the measured string tension a (from plane-plane correlators of Polyakov loops) 
with a calculated "theoretical" string tension using as input the measured monopole 
density p: 
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Here a^^ is given in accordance to eqs. (|TT],^2D via 
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Pi/3) 
MP) 



(29) 



and Pv is defined in eq. (|^). 

The circles shown take into account all "active" monopoles, i.e. isolated ones and those 
from charged monopole clusters which might be thought to be responsible for the string 
tension. The ratio is close to unity indicating the fact that the charged monopoles provide 
the major contribution to the string tension, as expected. Note that in the deconfinement 
phase the string tension is non-zero due to finite-size effects discussed below. The squares 
in Figure |^(a) are related to the ratio ( pSf ) in which all monopole are taken into account. 
In both phases this ratio is smaller than unity: a neutral fraction of the monopoles bound 
in the small dipole pairs does not contribute to the string tension. 

The small "string tension" remaining after passing the deconfinement transition at 
this finite lattice can be explained from the point of view of the dipole picture as follows. 
Test particles separated by distances smaller than sizes of certain dipoles are influenced 
by the constituent monopoles of those dipoles. The monopoles give contribution to the 
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string tension term. On the finite lattice the maximal distance between the test particles 
is of the order of the lattice size. Therefore, dipoles of the same size could be responsible 
for the non-vanishing small string tension. Dipoles of these sizes may really be present 
in the deconfinement phase (with a probability decreasing with the increase of the lattice 
size). This, however, does not contradict the criterion used to locate the phase transition 
in the previous Sections since the dipoles of such large size are heavily suppressed. 

The dipole formation due to Coulomb forces also happens at zero temperature. This 
effect increases the monopole density compared to that in the " bindingless" world. To 
check this we compare in Figure ^ the total density of the total density of monopoles and 
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Figure 8: The density of all monopoles and of monopoles in charged clusters vs. (3 for a 
32^ lattice compared to prediction (pi]). 



the exclusive density of monopoles residing in charged clusters (the latter includes free 
monopoles and anti-monopoles) for a 32^ lattice. The charged monopoles comprise around 
55% of the total monopole density. This ratio does not depend on the value of the coupling 
constant (3 indicating that the scaling behaviour of charged and neutral clusters is the 
same. The charged fraction of the monopoles is perfectly described by the "bindingless" 



formula (|TTD for the monopole density. This formula is incorporated implicitly into the 
theoretical prediction of the string tension (^2]) which works well according to Ref. ||T9| . 
Thus only the monopoles from the charged clusters (including separate monopoles) con- 
tribute to the string tension while the binding effect causes the appearance of a large 
fraction of "inactive" neutral clusters. 

Finally, we have measured the "spatial string tension" : the coefficient in front of the 
area term in the spatial Wilson loops. This string tension ag has been obtained by means 
of the standard diagonal Creutz ratios. The results are presented in Figure ^ as a function 
of p. As expected, the spatial string tension does not vanish and behaves smoothly across 
the deconfining phase transition. In contrast, we show in this Figure also the "true" 
temporal string tension extracted from temporal Wilson loops which drops down to the 
level of the finite-volume correction that we have just discussed. 

At sufficiently high temperatures the system might be treated as two-dimensional with 
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Figure 9: The spatial string tension vs. p. 



an effective 2D coupling constant P^'^^^ = Lt (3. Moreover, since the monopole density 
is low at large (3 (in deconfinement), the model becomes effectively non-compact. Thus 
the spatial Wilson loop behaviour in this regime is given by the perturbative one-photon 
exchange. In two dimensions the Coulomb law provides the linearly confining potential, 
V"(2^)(_R) = i?/2, corresponding to the spatial string tension. 
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(30) 



which is shown in Figure ^ by the dashed line. The spatial string tension data and the 
curve approach each other for sufficiently large (3. However, in the confinement phase the 
monopoles give a significant contribution to the spatial string tension. 



6 Summary 

In this paper we have considered a mechanism of the finite temperature deconfinement 
phase transition in three dimensional compact electrodynamics based on the monopole 
binding. The considerations are similar to those given in Ref. P] for the continuum 
theory and they incorporate features of the lattice geometry. This allows us to predict 
the pseudocritical coupling as a function of the lattice size. 

In our numerical simulations we have demonstrated that the monopoles are sensitive 
to the phase transition despite the fact that the monopole density itself behaves smoothly 
across the transition. The pseudocritical couplings found by the Binder cumulants of the 
density are very close to that identified using the Polyakov loop susceptibility. We stress 
that we did not intend to study the finite size scaling behavior of this model. 

Based on the observation to find (3c in this way we have studied the monopole properties 
in more detail. We have found that both the monopole density and the string tension differ 
from the predictions based on a model which does not take into account the monopole 
binding effects. However we have found numerically that the ratio between these two 
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quantities derived in that model (given by eq. (^)) remains valid in the confinement 
phase. 

We have observed that the dipole formation happens both in the confinement and 
deconfinement phases. In the deconfinement phase tightly bound dipoles — which are 
safely identified by a cluster algorithm — dominate in the vacuum. The dipoles are 
oriented dominantly in the temporal direction. These features are in agreement with 
general expectations discussed in the Introduction and in Section |^. 

At the confinement phase transition we observe mostly clusters with two constituents 
or single monopoles and anti-monopoles. Decreasing further the temperature (or (3), the 
monopoles become dense and form connected clusters (on a coarser and coarser lattice) 
inclosing various numbers of monopoles and antimonopoles. The largest clusters are more 
and more spherical. Whether the observed properties of the dipole gas formation survives 
in the continuum limit deserves an additional study. 

When the phase transition is mediated by charged objects, one could expect that 
external fields will infiuence the phase transition. In our case the natural question arises 
what will happen to the confinement-deconfinement phase transition. For non-Abelian 
theories in 3 + 1 dimensions it was recently concluded, from a study of the expectation 
value of the Polyakov loop that confinement seems to become restored under the 
infiuence of an external chromomagnetic field. In an accompanying paper |PT[ we will 



report on a study of our model under such external conditions, concerning the infiuence 
of confinement and relevant properties of the monopole system. 
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